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Abstract. We discuss the various mechanisms potentially at work in hadroproduction of heavy quarkonia 
in the light of computations of higher-order QCD corrections both in the Colour-Singlet (CS) and Colour- 
Octet (CO) channels and the inclusion of the contribution arising from the s-channel cut in the CS channel. 
We also discuss new observables meant to better discriminate between these different mechanisms. 



1 Introduction 

Heavy quarkonia are among the most studied bound quark 
systems. It is however fair to say that, for the time being, 
there is no consensus concerning which mechanisms are ef- 
fectively at work in their production in pp collisions, that 
is at the Tevatron and RHIC. By extension, available theo- 
retical predictions for the production rates at the LHC are 
rather uncertain. For recent reviews, the reader is guided 
to [TJ[2J[2] along with some perspectives for the LHC [I] . 

To what concerns the T states, the latest NLO pre- 
dictions [5] in the QCD-based approach of the Colour- 
Singlet Model (CSM) [B] including for the first time some 
of the important NNLO a| corrections show a satisfactory 
agreement with the data coming from the Tevatron [318]. 

On the other hand, to what concerns the charmonium 
family, it is well-known that the first measurements by the 
CDF Collaboration of the direct production of J/ip and ip' 
at s/s — 1.8 TeV [9"1[TU] brought to light a striking puzzle. 
They indeed found much larger rates than the prediction 
of the CSM. 

Since then, other approaches were proposed (e.g. the 
Colour-Octet Mechanism (COM) from NRQCD QT]) or 
revived (e.g. the Colour-Evaporation Model (CEM) [12] ). 
Those are unfortunately not able to reproduce in a con- 
sistent way experimental studies of both cross-section and 
polarisation measurements for charmonia at the Tevatron 
[U3"1[L4"] along with the cross sections measured by PHENIX 
at RHIC [15] . For instance, the seemingly solid prediction 
of the COM for a transverse polarisation of ip's produced 
at high transverse momentum is clearly challenged by the 
experimental measurements. The most natural interpreta- 
tion of such a failure of NRQCD is that the charmonium 
system is too light for relativistic effects to be neglected 
and that the quark- velocity expansion (v) of NRQCD [11 
may not be applicable for the rather "light" cc system. 
This might be indeed the case in view of the aforemen- 
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tioncd agreement between theory and the available exper- 
imental data on production in pp (and inclusive decays) 
of the significantly heavier T. In this case, relativistic cor- 
rections are expected to be less important and the leading 
state in the Fock expansion, i.e. the heavy-quark pair in a 
colour singlet 3 Si, to be dominant. This would in turn ex- 
plain why a computation incorporating the sole CS chan- 
nel (LO contribution in the w-expansion of NRQCD) is suf- 
ficient - when P^f contributions are considered though. 

Recently, many new theoretical results became avail- 
able. Some completed our knowledge of quarkonium pro- 
duction like the up-to-date proof [TB] of NRQCD factori- 
sation holding true at any order in v in the gluon-fragmen- 
tation channel - therefore relevant at large Pt where this 
channel should dominate and the QCD corrections which 
we shall discuss in the next section. Others mainly ques- 
tioned our understanding of the mechanisms at work in 
heavy-quarkonium production. Firstly, a complete survey 
of fixed-target measurements [17] showed that the univer- 
sality of the Long Distance Matrix Elements (LDMEs) 
of NRQCEQ cannot certainly be extended to the descrip- 
tion of low-Pp data. Secondly, NRQCD factorisation was 
shown to require modifications in fragmentation regions 
where 3 heavy quarks have similar momenta [181. Thirdly, 
the c-quark fragmentation approximation was shown [19 
to be only valid at much higher Pt than expected from 
the pioneering works |20j . 

Parallel to that, we investigated in [2X1122] on the cut 
contributions due to off-shell and non-static quarks. In 
particular, we questioned the assumption of the CSM^] 
that takes the heavy quarks forming the quarkonium (Q) 
as being on-shell [6]. If they are not, the usual s-channel 
cut contributes to the imaginary part of the amplitude 



Let us remember that the universality of the LDMEs 
is clearly challenged by HERA data on photoproduction of 
J/i> 0. 

2 Note on the way that this assumption is also present in 
NRQCD. 
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and need to be considered on the same footing as the 
CSM cut. A first evaluation [2T] of the latter incorpo- 
rating constraints for the low- and large- Pt (the scaling 
limit) region give rates significantly larger than the usual 
CSM cut. Moreover, low-Pr data from RHIC are very well 
described without need of re-summing initial-gluon contri- 
butions. However, as expected [2T], this approach under- 
estimates the cross-section at large values of Pt and other 
mechanisms have to be considered in this region. 

In section 2, we present the latest results available 
on QCD corrections to hadroproduction of J/tp, ip' and 
T(nS). In section 3, we discuss how the s-channel cut 
contribution to the CS channel can be evaluated and we 
present a comparison with data. In section 4, we briefly 
review other recent theoretical results. In section 5, we 
show how the study of the production of quarkonia in as- 
sociation with a heavy-quark pair of the same flavour may 
be used to disentangle between the different mechanisms 
proposed to explain quarkonium production. Finally, we 
present our conclusions and outlooks. 



2 QCD corrections 

More than ten years ago now, the very first NLO calcula- 
tion on quarkonium production to date became available. 
It was centred on unpolarised photoproduction of ip [23] 
via a colour-singlet (CS) transition. Later on, NLO cor- 
rections were computed for direct 77 collisions 24,25 for 
which it had been previously shown [25] that the LO CS 
contribution alone was not able to correctly reproduce the 
measured rates by DELPHI [27]. NLO corrections have 
also recently been computed for the integrated cross sec- 
tion of two J/V'-production observables at the P-factories: 
J/tp + cc [28] and J/tp + rfr- [25]. As of today, only the 
full colour-octet (CO) contributions to direct 77 collisions 
have been evaluated at NLO for P T > [241125] , 

At the LHC and the Tevatron, ip and T production 
proceeds most uniquely via gluon- fusion processes. The 
corresponding cross section at NLO (a| for hadroproduc- 
tion processes) are significantly more complicated to com- 
pute than the former ones and became only available one 
year ago 30, 19 . We shall discuss them in the next section. 

The common feature of all these calculations is the sig- 
nificant size of the NLO corrections, in particular for large 
transverse momenta Pt of the quarkonia for the computa- 
tions of differential cross sections in Pt- In jp an pp colli- 
sions, QCD corrections to the CS production indeed open 
new channels with a different behaviour in Pt which raise 
substantially the cross section in the large- Pt region. 

Let us discuss this shortly for the gluon-fusion pro- 
cesses which dominate the yield in pp. If we only take 
into account the CS transition to 3 Si quarkonia, it is well 
known that the differential cross section at LO as a func- 
tion of Pt scale like P^ 8 [BJ. This is expected from con- 
tributions coming from the typical "box" graphs of Fig. [T] 
(a). At NLO [3UlfT9] . we can distinguish three noticeable 
classes of contributions. First, we have the loop contribu- 



tions as shown on Fig. [l] (b), which are UV divergent 3 
but as far their Pt scaling is concerned, they would still 
scale like Pf 8 . Then we have the i-channel gluon ex- 
change graphs like on Fig. [I] (c). They scale like P^ 6 . 
For sufficiently large Pt, their smoother Pt behaviour 
can easily compensate their as suppression compared to 
the LO (a|) contributions. They are therefore expected 
to dominate over the whole set of diagrams up to a s . 
To be complete, we should not forget the a s contribu- 
tions from Q+QQ (where Q is of the same flavour as the 
quarks in Q). Indeed, one subset of graphs for Q + QQ 
is fragmentation-like (see Fig. [I](d)) and scales like P^ 4 . 
Such contributions are therefore expected to dominate at 
large Pp, where the smoother decrease in Pt is enough 
to compensate the suppression in as and the one due to 
the production of 4 heavy quarks. As mentioned above, in 
practice [IS], this happens at larger P T than as expected 
before [2U]. We shall come back to this channel later. In 
the next sections, we shall discuss the impact of the NLO 
corrections to the CS channels and then a first computa- 
tion including the a priori dominant a s contributions i.e. 
topologies illustrated by Fig. [T](e) and (f). 




(d) (e) (f) 




Fig. 1. Representative diagrams contributing to 3 Si hadropro- 
duction via Colour-Singlet channels at orders a| (a), a| 
(b,c,d), a| (e,f) and via Colour-Octet channels at orders a| 
(g,h). The quark and antiquark attached to the ellipsis are 
taken as on-shell and their relative velocity v is set to zero. 

To what concerns the CO contributions, the effects of 
NLO (here a s ) contributions are expected to be milder. 

3 These divergences can be treated as usual using dimen- 
sional regularisation, see e.g. [30] . 
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Indeed, the contributions from CO transitions are a pri- 
ori suppressed by v 4 and the reason why they can still 
appear significant comes from a lower power in as for 
similar topologies (and thus similar Pt scaling). For in- 
stance, compare Fig. [I] (g) to Fig. [I] (e) and Fig. [I] (h) 
to Fig. [j] (c). At LO, Pf 6 and Pf 4 scaling arc therefore 
already present. As a result, including ctg contributions 
will not open any new channels and NLO corrections are 
expected to be described by a roughly constant K factor. 
As we shall see, this is indeed the trend seen in the results 
of [3T] in which the NLO corrections to l Sf ] and 3 Sj 81 
colour octets going to J/tp were considered. 



2.1 NLO corrections for Colour-Singlet channels 

Let us first present a comparison between the measure- 
ments by the CDF collaboration and the result for the 
J/tp obtained following the procedure explained in |30l 
119] . It is worth noting that the Xc cross sections are not 
available for now at NLO accuracy. This would be neces- 
sary if we wanted to predict at this accuracy prompt- J/tp 
production cross section, in order to compare with the 
most recent measurements of RUN II [32] . These focused 
only on the prompt yield. As a makeshift, we have mul- 
tiplied those data by the averaged fraction of direct J/ijj 
measured during RUN I |10j for the rather similar beam 
energy 1.8 TeV and a similar range in Pt and rapidity: 
(F dlrect ) = 64 ± 6%. 

In our calculation, we set m c = 1.5 ± 0.1 GeV and 
m b = 4.75 ± 0.25 GeV. We used the PDF set CTEQ6L1 
(resp. CTEQ6JVI) [33J for LO (resp. NLO) cross sections, 
and always kept the factorisation scale equal to the renor- 
malisation scale: (if = fi r . Except for the associated pro- 
duction channel where we took /io = J (2m,Q) 2 + P T , the 

central scale is fixed at //o = \J + Pp and then was var- 
ied by a factor of 2. To what concerns the non-perturbative 
inputs, we used the values related to the BT potential [33] : 

(O j /*?S 1 i ] )) = 1.16 GeV 3 and (0 T ( 3 s[ 1] )) = 9.28 GeV 3 . 
For Y(1S) production, we considered the prompt measure- 
ment at s/s =1.8 TeV in [7 , multiplied by the averaged 
direct fraction obtained in [35]: (F direct ) = 50 ± 12%. 

In both case, we have an illustration of the previous 
discussion. The differential cross section for the LO con- 
tribution, i.e. gg — > J/tpg, has the steepest slope and is 
already an order of magnitude smaller than the NLO con- 
tribution at Pt — 10 GeV. The differential cross section 
for Q + QQ has the smoothest slope. In the case of J/tp, 
it starts to be significant for Pt > 20 GeV. For the T, the 
suppression due to the production of 4 b quark is stronger 
and this yield remains negligible in the accessible value of 
Pt- The bands denoted NLO refer to all the contributions 
up to order ctg. 

Those results were recently confirmed in [36, 37] . In the 
latter papers, the polarisation information was kept and 
the observable a was also computed. However, it is impor- 
tant to stress that for tp and T production the CS yields 
predicted at the NLO accuracy are still clearly below the 
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Fig. 2. Differential cross sections at NLO accuracy as function 
of the quarkonium transverse momentum Pt at the Tevatron 
(V5= 1-96 TeV). 



experimental data especially at large Pt- In this respect, 
the predictions for the polarisation at this order cannot 
be usefully compared to the data. 

The conclusion is that, in general, the inclusion of NLO 
contributions bring the CS predictions considerably closer 
to the data, although agreement is only reached at NLO 
in the photoproduction case [23] . 



2.2 NLO corrections for Colour-Octet channels 

As aforementioned, NRQCD has reached a certain suc- 
cess by explaining the main features of charmonium and 
bottomonium hadroproduction via the introduction of the 
Colour-Octet (CO) mechanism. It indeed provides a good 
description of the Pr-differential cross-section for the di- 
rect J/tp and tp' for Py > 5 GeV as measured by CDF 
in pp pntj], A reasonable agreement was also obtained 
with the first PHENIX measurements in pp at y/s = 200 
GeV [351l39j . In both cases, the cross-section is dominated 
by the gluon fragmentation into a colour-octet S'-wave 
state. Following the heavy-quark spin symmetry |11] of 
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Fig. 3. Full computation at NLO for (Left) Y(IS) + X (dashed line) vs. Y(1S) + 2 light partons with a cut on s™ ln (grey 
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NRQCD, the latter mechanism leads to transversally po- 
larised J/ip and ip' , the parent fragmentating gluon being 
mostly on-shcll and thus transversally polarised at high 
Ft- 

However, J ftp and ip 1 are not seen to be transverse by 
the CDF experiment [IS]- It measured a slight longitudi- 
nal polarisation for both the prompt J/ip and direct tp' 
yield. It is worth noting here that the feed-down from \c 
can influence significantly the polarisation of the prompt 
J/ip yield - this was taken into account in the NRQCD- 
based predictions [40 . Moreover, the recent preliminary 
result from PHENIX (41] indicates a polarisation compat- 
ible with zero for the total J/ip production at forward 
rapidity (1.2 < \y\ < 2.2), but with large uncertainties. 

Very recently, CO contributions from S waves ( 1 ^q 8 ' 

and 3 sj 8 ^) have become available [31] for hadroproduc- 
tion. A complete phenomenological study is not yet avail- 
able though. Anyhow this confirms that NLO corrections 
do not affect significantly the Pt dependence as expected 
from the introductory discussion of this section. 

Let us define K factors as the ratios of NLO to LO 
cross section for a given CO channel. For the Tevatron, 
they are about 1.2 for the state and 1.1 for the 3 sj 8 ' 
(at the LHC, they are both about 0.8). Consequently, the 
value of the CO Long Distance Matrix Elements (LDMEs) 
fit to the Tevatron data at LO (O j ^ ( 3 sf ] )) ~ 0.0012 

GeV 3 and (0 J ^(/Sf ] )) ~ 0.0045 GeV 3 would be at 
most reduced by 15%. In this respect, the NLO corrections 
to the octets do not improve the universality of the ma- 
trix elements when the idea of the dominance of the CO 
transitions is confronted to the data on photoproduction 
from HERA. 

According to the author of [31], it is not possible to 
obtain a satisfactory Pt distribution in terms of a unique 
(On) value when considering the whole range in Pt anal- 
ysed by CDF. More precisely, they did not consider the 
experimental data with Pt < 6 GeV, for which it seems 
that other mechanisms have to be at work if we believe 



that the COM is responsible for the major part of the 
cross section at large Pt- 

This in any case emphasises the need for more work 
dedicated to the description of the low-Py region and 
maybe the relevance of the study of s-channel cut con- 
tribution, which we discuss later. Last but not least, the 
polarisation from CO transitions appears not to be modi- 
fied at NLO with respect to LO results. Overall, this recent 
first study of CO contributions at NLO in hadroproduc- 
tion at Pt > sounds like a confirmation of the flagrant 
discrepancy between the NRQCD predictions for the po- 
larisation of the J/ip and the experimental measurements 
from the CDF collaboration [13 . 

2.3 QCD corrections up to a| 

As noted above, the discrepancy between the NLO com- 
putations for the CSM and the experimental data, both 
for ip and T still grows with Pt . If we parallel that to the 
existence of new P^ 4 channel at order a|, it is reasonable 
to wonder what their size are effectively . 

In fact, their contributions can be evaluated in a rel- 
atively "simple"]^] and reliable way by computing the a| 
contributions consisting in the production of a Q with 3 
light partons (noted j thereafter). Among them are the 
topologies of Fig. [T] (d) (gluon fragmentation) and Fig. [I] 
(e) ("high-energy enhanced"), these close the list of kinc- 
matical enhancements from higher-order QCD corrections. 
This a| subset being the LO for a physical process (pp — > 
Q + jjj), its contribution is finite except for soft and 
collinear divergences. 

To avoid such divergences, we impose a lower bound on 
the invariant- mass of any light partons (s.y). For the new 
channels opening up at a|, and which specifically interest 
us, the dependence on this cut is to get smaller for large 
Pt since no collinear or soft divergences can appear there. 

4 "simple" compared to a full -out-of-reach- NNLO compu- 
tation and thanks to the automated generator of matrix ele- 
ments MadOnia [42]. 
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For other channels, whose LO contribution is at ctg or a.g, 
the cut would produce logarithms of s^/s™ 11 . Those can 
be large. Nevertheless, they can be factorised over their 
corresponding LO contribution, which scales at most as 
P^ 6 . The sensitivity on s™ m is thus expected to come to 
nothing at large Pt- 

Thanks to the exact NLO computation of [3D], such 
a procedure can be tested for the process pp — ► Q + jj. 
For instance, the differential cross section for the real at 
corrections, T(1S) + jj production, is displayed in Fig. [3] 
(Left). The grey band illustrates the sensitivity to the 
invariant-mass cut s™ ln between any pairs of light partons 
when it is varied from 0.5m^ to 2m%. The yield becomes 
insensitive to the value of s™ m as Pt increases, and it re- 
produces very accurately the differential cross section at 
NLO accuracy. In the charmonium case, the similar con- 
tributions from pp — > ip 1 + jj matches even better, for 
lower Pt and with a smaller dependence of s™ the full 
NLO computation, as seen on Fig. [3] (Right). 



We now turn to the results concerning the real con- 
tribution at do, which we refer to as NNLO*. We used 
the approach described in Ref. [32], which allows the au- 
tomatic generation of both the subprocesses and the cor- 
responding scattering amplitudes. The differential cross- 
sections for Y(1S) and ip(2S) are shown in Fig.Pl] The red 
band (referred to as NNLO*) corresponds to the sum of 
the NLO yield and the Q + jjj contributions. In the T 
case, the contribution from T with three light partons fills 
the gap between the data and the NLO calculation, while 
for the Tp(2S) there seems to remain a small gap between 
the NNLO* band and the preliminary CDF data [43]. In 
both cases, the a| contribution is very sensitive to the 
choice of the renormalisation scale, [i r . This is expected: 
for moderate values of the Pt, the missing virtual part 
might be important, whereas at large Pt, the yield is dom- 
inated by Born-level a|-channels from which we expect a 
large dependence on \i r . Even though the uncertainty on 
the normalisation is rather large, the prediction of the Pt 
shape is quite stable and agrees well with the behaviour 
found in the data [71I51FB], 
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Fig. 4. Comparison between differential cross sections at NLO 
and NNLO* accuracy as function as function of the Q trans- 
verse momentum Pt at the Tevatron (^/s = 1.96 TeV) and the 
data for (a) T{1S) [7] & (b) direct ip(2S) g3]. 
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Fig. 5. Polarisation of (a) T(nS) ((b) ip(2S)) directly produced 
as function of its transverse momentum Pt at the Tevatron. 



Concerning the polarisation, the direct yield is pre- 
dicted to be mostly longitudinal, see Fig. [5] (a). However, 
existing experimental data for T are centred on the prompt 
yield [71l44j . In order to draw further conclusions, we would 
need first to gain some insights on NLO corrections to P- 
wave production at Pt > . Yet, since the yield from 
P-wave feed-down is likely to give transversely polarised 
T, the trend is more than encouraging. To what concerns 
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ip{2S), one should be very careful before any comparison 
with experimental measurements since the yield is not ex- 
actly reproduced. Having this in mind, one sees in Fig. [5] 
(b) that the trend for longitudinally polarised tp(2S) is 
reproduced but more marked. At very large Pt where the 
contribution from tp(2S) + cc becomes more and more sig- 
nificant the polarisation gets slightly less negative. In any 
case, further investigations are needed to draw any con- 
clusions. 



3 s-channel cut contribution 

In this section, we briefly review our first evaluation of the 
s-channel cut contribution in hadroproduction of J/ip |21l 
122] and present outlooks for necessary future investiga- 
tions in this direction. 

If the quarks which constitute the Q are not on their 
mass shell, it is not possible to factorise in a gauge invari- 
ant way the amplitude for the production of those quarks 
and the one responsible for their binding into the Q. In 
other words, if the latter amplitude is given by a (3-point) 
Bcthe-Salpeter vertex function, the set of Feynman dia- 
grams shown on Fig. [6] (a,b) constructed from this vertex 
will not be gauge-invariant. Indeed, we have to introduce 
a 4-point function (or contact term) and build up from 
it new contributions to the production amplitude (Fig. [6] 
(c)). Note that the 4-point function ccJ/ipg could be in- 
terpreted as a dynamical generalisation of a CO matrix 
elements for the transition between a C — +1 CO state 
and a J/ip. 




Fig. 6. (a) & (b) Leading-order (LO) s-channel cut diagrams 
contributing to gg — > Qg with direct and crossed box diagrams 
employing the 3-point ccQ vertex. The crosses indicate that 
the quarks are on-shell. (c) Box diagram with the (4-point) 
ccQg contact term mandated by gauge invariance. 



where g s is the strong coupling constant, Tf k the colour 
matrix, and // and v are the Lorentz indices of the outgoing 
J /if) and gluon, respectively. For simplicity, we have sup- 
pressed all indices on the left-hand side. The ccJ/ip vertex 
function p( 3 ) with the kinematics of the direct graph is 
denoted here by Zi and for the crossed graph by P 2 , 
P x = P (ci - f , P) and P 2 = P (c 2 + f , P) , as shown in 
Figs. |6ja) and (b). 

One easily verifies that any contact current defined as 
followf0 



MX 



(2c 2 + qr(r 1 -F) (2ci - q) v (P a - F) 



( c 2 + q) 2 — m 2 



(ci - qf 



(3) 



satisfies the gauge-invariance requirement |21j for any value 
of the function F — P(c 1; 02,(7) ■ 

Nonetheless, the function F[c.\,C2,q) must be chosen 
so that the current ^ satisfies crossing symmetry (i.e., 
symmetry under the exchange c\ <-> — c 2 ) and is free of sin- 
gularities. The latter constraint implies F = Po at either 
pole position, i.e., when (c 2 + <7) 2 = m 2 or (ci — q) 2 = m 2 , 
where the constant Po is the (unphysical) value of the 
momentum distribution P(p, P) when all three legs of the 
vertex are on their respective mass shells. In principle, em- 
ploying gauge invariance as the only constraint, we may 
take F = Po everywhere. This corresponds to the min- 
imal substitution discussed by Drell and Lee [13] (for a 
complete derivation see |50j ) who pointed out, however, 
that this does not provide the correct scaling properties 
at large energies, which means within the present context 
that F — Po would not lead to the expected Pt scaling of 
the amplitude. See |51j for a numerical comparison with 
data. 

In order to obtain a correct scaling at large Pt and a 
behaviour close to the minimal substitution at low Pt, we 
have chosen [2T1I5T] 



F(c 1 ,c 2 ,q) = P - h{ci ■ c 2 ) 



(Po - PQ (Pp - P 2 ) 
Po 



(4) 



where the (crossing-symmetric) function h(c\ ■ c 2 ) rises 
to become unity for large relative momentum. The phe- 
nomenological choice for the interpolating function h used 
in our calculations is 



Gauge invariance relates the 4-point function to the 3- 
point one, but not univocally since it does not constrain its 
transverse contribution with respect to the emitted gluon. 
However there is an elegant way to parametrise this free- 
dom by an auxiliary function P |45ll46 47.48], To see this, 
let us define the 3-point function 



P( 3 >(P,P) = P(P,P) 7M 



(1) 



where P = p\ — p 2 and p = (jp\ + p 2 )/2 are the total and 
relative momenta, respectively, of the two quarks bound as 
a quarkonium state, with pi and p 2 being their individual 
four-momenta and the 4-point one by 



P< 4 > - -ig s T? k M» cl » 



(2) 



h{c\ ■ c 2 ) = 1 — a- 



(ci • c 2 + m 2 ) 



(5) 



with two parameters, a and k. 

Figure [TJa) shows our results for y/s — 1.8 TeV in 
the pseudorapidity range \rj\ < 0.6 with parameter values 
a = 4 and k — 4.5 GeV fixed to reproduce, up to P T — 
10 GeV, the cross-section measurement of direct J/ip by 
CDF [TO], the usual LO CSM from gg -> J/ipg |] and 
LO CSM from gg -> J/ipcc [T5]. Our results fit well the 
CDF data up to about P T = 10 GeV. At higher P T , our 
curve falls below the data as expected from the genuine 



5 We have taken c\ = c% = m 2 and P 2 = M 2 with m and M 
being the masses of the quark and the J/tp. 
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Fig. 7. (a) Comparison between polarised (or and ah) and unpolarised (a to t) cross sections [with parameters a = 4, k — 4.5 
GeV in Eq. (JsJ ] , LO CSM contributions, and CDF experimental data |10j at the Tevatron (^/s — 1.8 TeV, pseudorapidity 
|r/| < 0.6); (b) Comparison between a T , a L , a to t and PHENIX data [15] at RHIC (s/s = 200 GeV, rapidity \y\ < 0.35). Taken 
from [21]. 



l/P T scaling of a LO box diagram (see the discussions of 
the previous section). 

It is interesting to note the different Pt behaviours of 
<jt and <jl leading to a dominance of the latter at large Pt 
and a negative value for the polarisation a |21j at mid and 
large Pt- Figure [TJb) shows our results at %fs = 200 GeV, 
still with a — 4 and k — 4.5 GeV, compared with the 
PHENIX data [J5]. 

Through this first evaluation of the s-channcl cut con- 
tribution to the imaginary part of the production ampli- 
tude, incorporating low- and large-energy constraints as 
well as gauge invariance, we have shown that this cut can 
be significant. It is even possible to obtain a very good 
fit of the data from CDF at mid Pt by proper choices 
of the parametres of our 4-point function. With the same 
parametres, we obtained an excellent description of the 
data taken at RHIC and down to very low Pt without re- 
summing initial-state gluon contributions. The s-channcl 
cut indeed has a threshold at low s (thus low Pt) which 
corresponds to the energy needed to put the two c-quark 
on-shell. 

Now that we have seen that the s-channel cut mat- 
ters at low- and mid-Py, it is necessary to have in the 
future a first evaluation of the contribution of the real 
part itself. On the other hand, we can start testing our 
parametrisation of the 4-point function, in photoproduc- 
tion for instance, or in any other process involving a final 
state gluon. 



4 Other theoretical advances 

Beside the theoretical advances concerning QCD correc- 
tions and the inclusion of the s-channel cut contribution 
discussed in the previous sections, several interesting the- 
oretical results have been obtained in the recent years. Let 
us review some of the most significant ones briefly. 

On the side of NRQCD, Nayak, Qiu and Sterman pro- 
vided an up-to-date proof [IB] of NRQCD factorisation 
holding true at any order in v in the gluon-fragmentation 



channel. They showed that new definitions of NRQCD 
matrix elements incorporating QCD Wilson lines were to 
be used, but that this was not to affect the existing phe- 
nomenological studies. 

Last year, Collins and Qiu |52j showed that in gen- 
eral the fey-factorisation theorem does not hold in pro- 
duction of high-transverse-momentum particles in hadron- 
collision processes, and therefore also for ip and T. This is 
unfortunate since many studies [53, 54 . 55 , 56 , 57 , 58 , 59 , 60 , 
|6"T1I6"2"] . predicting mostly longitudinal yields and smaller 
CO LDMEs, in better agreement with the idea of LDME 
universality, were based on the hypothesis of such a fac- 
torisation in hadroproduction. 

Besides, the c- and ^-fragmentation approximation was 
shown to fail for the Pt ranges accessible in experiments 
for quarkonium hadroproduction. By studying the entire 
set of diagrams contributing to ip and T production in 
association with a heavy-quark pair of the same flavour, 
we have shown [TIT] that the full contribution was signif- 
icantly above (typically of a factor of 3) that obtained 
in the fragmentation approximation. A precision of 10% 
accuracy, say, can only be obtained at very large Pt ■ 
P T > 60 GeV for V and P T > 100 GeV for T. Note that 
the same observation was previously made for the process 
77 — > J/ipcc [53] and also for the B* hadroproduction, for 
which it was noticed that the fragmentation approxima- 
tion was not reliable at the Tevatron [M]l65]. 

Moreover, still in double-heavy-quark-pair production, 
the notion of colour-transfer enhancement was introduced 
by Nayak, Qiu and Sterman [18] . If three out of the four 
heavy quarks are produced with similar velocities, then 
there is the possibility that colour exchanges within this 3- 
quark system could turn CO configurations into CS ones, 
thus could effectively increase the rate of production of 
CS pairs. They finally discussed the introduction of spe- 
cific new 3-quark operators -beyond the usual ones of 
NRQCD- necessary to deal with such an issue. A study 
of the colour-transfer effects in hadroproduction is still 
awaited for. 
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5 Associated production channels 

As previously discussed, the results of QCD corrections 
for Y production seem to indicate that the CS transitions 
are dominant. Eventually, this should put an end to the 
controversy related to T production. Contrariwise, the sit- 
uation remains unclear for the charmonium case. NLO 
corrections complemented by some dominant ct| contri- 
butions are large and seem to bring the prediction for 
the CS transitions very close to the data in the ip' case 
for instance (see Fig. H (b)). Yet, theoretical uncertain- 
ties remain large and there seems to be some space left 
for CO transitions. Careful comparisons are still therefore 
due with polarisation observables. In this case, the theo- 
retical uncertainties would certainly be competitive with 
experimental ones, for instance on prompt J/ip yield [66 . 
This requires however some knowledge on the QCD correc- 
tions to the P-wave CS yield. For the time being, nothing 
is known on this side. 

It is therefore vital in order to progress in the under- 
standing of the mechanisms responsible for heavy quarko- 
nium production to introduce, compute and measure new 
observables. One of those is the hadronic activity around 
the quarkonium 67J. Historically, UA1 compared their 
charged-track distributions with Monte Carlo simulations 
for a J/ip coming from a B and a, J ftp coming from a 
Xc [68ll69j. At that time \c feed-down was still expected 
to be the major source of prompt J /ip- Following either 
the idea of CO transitions or of CS transitions at higher- 
orders, we however expect now more complex distributions 
even for the prompt yield. It is therefore not clear if such 
methods are suitable to size up the P-feeddown otherwise 
than with the measurements of a displaced vertex typical 
of a B decay. 

We therefore urgently need observables rather easy to 
predict and likely to test the many production models 
available [T|[2] . We argue here that the study of associated 
production channels, first in pp collisions, then in pA and 
AA, fills both these requirements. By associated produc- 
tion channels, we refer to ip + cc and T + bb. 

A further motivation for such studies is that simi- 
lar studies carried at P-factories showed an amazingly 
large fraction of J/ip produced in association with an- 
other cc pair. Indeed, the Belle collaboration first found 

a[ t+ e e-Z J j'/tX C x) t0 be 0.59+g;li±0.12. Thereafter, the 



Effl 

analysis was improved and they obtained |71j 



er(e + e -^J/i]j + cc) 
cr(e+e- ~* Jjf+X) 



0.82 ±0.15 ±0.14, 



(0) 



> 0.48 at 95% CL. 



Whether or not such a high fraction holds for hadropro- 
duction as well, is a question which remains unanswered. 
Analyses at the Tevatron (CDF and D$) and at RHIC 
(PHENIX and STAR) are already possible. As computed 
in PI] for the RUN2 at the Tevatron at ^fs = 1.96 TeV, 
the integrated cross-section are significant : 



a(J/ip + cc) x B(£ + £-) ~ 1 nb 
a{T + bb) x B(£ + £-) ~ 1 pb 



(7) 



As an illustration of the potentialities at RHIC, we 
chose to display in Fig. [8] the differential cross section 
for pp — > J/ip + cc computed for the STAR kinematics. 
Such studies could for instance be carried out by STAR 
in the next run with an integrated luminosities of around 
50 pb -1 if dedicated triggers are available [72] , 
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Fig. 8. Differential cross section for pp — + J/ip + cc as function 
of the J/ip transverse momentum Pt for the STAR kinematics 
(^= 200 GeV, \y\ < 1.0). 



Without taking into account the likely reduction of the 
CO LDMEs induced by the QCD corrections mentioned 
in the previous sections, the integrated cross sections were 
found in |73j to be dominated by the CS part, similarly to 
the differential cross section in Pt up to at least 5 GeV 
for ip and 10 GeV for T. In other words, such observables 
can be thought of as a test of the CS contribution, for 
the first time since the introduction of the idea that CO 
transitions would be the dominant mechanism responsible 
for quarkonium production at high transverse momentum. 
If the effect of CO transitions is confirmed to be negligible 
for the T, the T produced in association with a bb pair are 
predicted to be strictly unpolarised, for any Pt (see Fig. |9| 
for the LHC. 

Beside the property of discriminating between the CO 
and the CS transitions, the yield of ip in association with cc 
should show an a priori completely different sensitivity to 
the Xc feed-down than the inclusive yield. The same holds 
for T with bb with the Xb feed-down. To what concerns 
CS transitions, the P-wave yield is expected to be smaller 
than the S'-wave one, since they are being suppressed by 
powers of the relative velocity v and here there is no extra 
gluon needed to be attached to the heavy-quark loop to 
produce ip (or T) compared to P-waves as it is the case in 
the inclusive case. 

For the CO transitions, associated production Xc + cc 
can occur via the process gg — > gg for which the two final- 
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Fig. 9. Polarisation of an T produced in association with a bb 
pair at the LHC for y/s = 14 TeV for \y\ < 0.5. 



state gluons split into a cc pair, one of them hadronising 
into a Xc via the CO mechanism. This contribution is cer- 
tainly suppressed up to P T ~ 20 GeV. For larger P T , a 
dedicated calculation is needed. However, this mechanism 
would be very easily disentangled from the CS contribu- 
tions since both c quarks are necessarily emitted back to 
back to the \c and thus to the J/ip. 

Concerning the non-prompt signal, it would originate 
as usual from gg — > bb, where one b quark hadronises 
in if>. Usually, this hadronisation of the b produces the 
tf> with light quarks only. This means than we have one 
single c quark in the event. It is produced from the decay 
of the recoiling b quark and is therefore back to back to 
the if). The non- prompt signal would then be simply cut 
down by searching for a D meson near the if). Now it can 
happen that the hadronisation of the b produces the if) 
and a D meson. In this case, kincmatical cuts would not 
help to suppress the non-prompt yield. Fortunately, this 
is a priori suppressed compared to the first case and even 
more than the direct yield since there is here no gain in the 
Pt dependence since both gg — > bb and gg — > tp + cc cross 
sections scale like P^f 4 . A cross check by sizing up the non- 
prompt yield with a displaced vertex measurement would 
be anyhow surely instructive. 

Let us also mention that associated production has 
also been studied in direct 77 collisions in Ultra-peripheral 
collision (UPC) [74,. At least for direct 77 collisions, as- 
sociated production is the dominant contribution to the 
inclusive rate for Pt > 2 GeV/c. 

To conclude, studies can be carried on by detecting 
cither the "near" or "away" heavy-quark with respect to 
the quarkonia. There are of course different way to detect 
the D, B, or a 6-jet, ranging from the use of a displaced 
vertex to the detection of their decay in e or fi. As dis- 
cussed above, this has to be considered by also taking 
into account the different backgrounds. The forthcoming 
Quarkonium-event-generator Madonia 2 |75j will surely 
be of a great help to achieve this task. In any case, we 
hope that such measurements would provide with clear 
information on the mechanisms at work in quarkonium 
production^] 



6 Conclusion 

Recently significant progresses have been made in the eval- 
uation of the QCD corrections to quarkonium production. 
The situation sounds now rather clear for the bottomonia 
where an agreement has been eventually obtained using 
only CS channels when dominant a| contributions are 
incorporated. The polarisation predictions for the latter 
cases seem also quite encouraging considering CDF [7j and 
D% [44] measurements. Yet, confirmations are awaited for 
from the LHC. 

On the other hand, those or s contributions could be 
still unable to bring agreement with the measured Pt- 
diffcrcntial cross-section of the direct charmonia. Ded- 
icated further studies are however needed especially to 
what concerns the feed-down from P-waves which is not 
known at NLO accuracy. In the charmonium case, we have 
also seen that s-channel cut can bring a significant con- 
tribution to the cross section at low Pt and hence a first 
evaluation of the real part of the production amplitude is 
needed. 

Additional tests are now undoubtedly needed beyond 
the mere measurements of inclusive cross section and po- 
larisation at the LHC. For instance, the hadroproduction 
of J/if) or T with a heavy-quark pair [19,73 appears to 
be a new valuable tool to separately probe the CS con- 
tribution, at least dominant at low-P^ (below 15 GeV), 
as well as the study of the hadronic activity around the 
quarkonium. 
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